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ABSTRACT 

The CMB maps obtained by observations always possess domains which have to be 
masked due to severe uncertainties with respect to the genuine CMB signal. Cosmo- 
logical analyses ideally use full CMB maps in order to get e. g. the angular power 
spectrum. There are attempts to reconstruct the masked regions at least at low reso- 
lutions, i. e. at large angular scales, before a further analysis follows. In this paper, the 
quality of the reconstruction is investigated for the ILC (7yr) map as well as for 1000 
CMB simulations of the ACDM concordance model. The latter allows an error esti- 
mation for the reconstruction algorithm which reveals some drawbacks. The analysis 
points to errors of the order of a significant fraction of the mean temperature fluc- 
tuation of the CMB. The temperature 2-point correlation function C($) is evaluated 
for different reconstructed sky maps which leads to the conclusion that it is safest to 
compute it on the cut-sky. 

Key words: Methods: data analysis, statistical; Cosmology: cosmic microwave back- 
ground, large-scale structure of Universe 
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1 INTRODUCTION 

The cosmic microwave background (CMB) provides one of 
the cornerstones of the cosmological concordance model. 
The statistical properties of our cosmological models have 
to match those of the CMB in order to give an admissible 
model. Thus, it is of utmost importance to reliably extract 
the statistical properties of the CMB. A main obstacle is the 
foreground emission of our galaxy and of other sources which 
restrict the area of the sky available for a sufficiently clean 
CMB signal, i. e. the full sky CMB signal has to be masked. 
One path of statistical analysis leads to the Fourier space 
in which the CMB is decomposed with respect to spherical 
harmonics Yj m (n) where the masked sky leads to a coupling 
between the Fourier modes since no full sky CMB is avail- 
able. In this paper, we do not delve into these difficulties, 
where upon an extensive literature exists, but instead follow 
the alternative path which allows an analysis directly in the 
pixel space, which in turn is more adapted to a masked sky. 
This analysis is based on the temperature two-point corre- 
lation function C(#), which is defined as 



C(<&) ~ (ST(h)ST(n)) with n-n = cos ■& 



(1) 



where 5T(h) is the temperature fluctuation in the direction 
of the unit vector h. The most direct way to deal with a 
mask is just to use only those pixels which are outside the 
mask. In this way it was discovered by the COBE team 
IjHinshaw et al.l 1 19961 ) that the correlation function C(#) 
possesses surprisingly low power at large angles i? > 60°. 



A surprising observation is made by using the ILC 
map, which represents a full sky CMB map obtained by 
the WMAP team (Gold et alj|20ich . Computing the corre- 
lation function C(#) using a mask leads to a correlation 
function having very low power at large scales, whereas 
using the complete ILC map leads to a correlation func- 
tion which possesses higher large sca le power b eing com- 
patible with the concordance model (|Spergel et all I2003T ) . 
One has to decide which correlation function C(#) corre- 
sponds to the true CMB sky: the one which is based on 
the safe pixels, i. e. those outside the mask, or the other 
one, which would imply that most of the large scale power 
is generated by those areas hidden by the galaxy, i. e. 
by those pixels which have experienced much larger cor- 
rections. This question has recently s timulated much dis- 
cussions, e.g. JCopi et all 120071. I2009J. l20ld : iHaiianl 120071; 
Aurich et alIl200Sl, l20ld: ISarkar et alj|2010l; iBennett et ail 
20ld ; lEfstathiou et al.ll2009l ; IPontzen and Peiri j|20ld ) . 



lEfstathiou et al.l |2009l ) emphasise that one has to start 
with the cut sky, but before the correlation function C(#) is 
computed the low-order multipoles have to be reconstructed. 
In this way, a stable result is obtained as long as the mask 
is not too large. The obtained correlation function is then 
the one with large power at large scales. For the detail s 
of the method see de Oliveira-Costa and Tegmarkl (2006); 
iBielewicz et al.l (J200J). Here we summarise only the most 
important ingredients. The data vector x containing only 
the pixel values outside the mask is related to the spherical 
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harmonic coefficients ai m represented as a by 
x — Y a + n , 



(a) KQ85 (7yr) mask, iV aidc = 512 



(2) 



where Yy denotes the corresponding values of Yi m . (hi) and 
h the noise. In the low-order multipole reconstruction, only 
the multipoles with I ^ Z max = 10 ... 20 are taken into ac- 
count. The methods of reconstruction differ by the choice of 
a square matrix A which determines the reconstructed a T 
by 



= (Y T AY)- 1 Y T Ax 



(3) 



Setting the matrix A equal to the unit matrix leads to 
the method of "direct inversion". To take the correla- 
tions between the pixels into account, one can choose 
the covariance matrix A -1 —< x ■ x T > for a recon- 
struction up to l max , which leads to the method used by 
Ide Oliveira-Costa and Tegmarkl |2006l ) with 



a: 



E 



21 + 1 
4ir 



Pi(hi 



)C, 



(4) 



ignoring the noise contribution. The sum in Q runs only 
over those multipole moments Ci, I > Z max , that are not to 
be reconstructed. 

The reconstructed sk y map up to l ma x i s then obtained 
from the coefficients a r . Ef stathiou et al.l (2009) compute 
the correlation function C(#) from such reconstructed maps 
for £ max = 5 to Z max = 20 by using the KQ85 or the KQ75 
mask provided by the WMAP team, and it is found that 
these C(i?) only show a negligible variation. This is inter- 
preted as a sign for a stable method to compute the tem- 
perature correlation in the presence of masks. It should be 
emphasised that the WMAP team bas es its 7 year investi - 
gation with respec t to CMB anomalies ([Bennett et al.ll2010P ) 
on the method of (Efstat hiou et al.ll2009l ). This leads to the 
conclusion that the correlation function agrees well with the 
ACDM concordance model and displays thus no anomalous 
behaviour. 

The system of equations © is over-determined, since 
there are much more pixel values than multipole moments 
with / ^ Z max as long as the mask and i max are not too large. 
Thus, there is the hope that the pixels outside the mask al- 
ready determine t he low-order multipoles . The caveat in the 
demonstration of (jEfstathiou et al.l 120091 ) is, however, that 
a large Gaussian smoothing of 10° is applied to the maps 
before the cut is applied. Due to this smoothing there is a 
transfer of information from pixel values within the mask to 
those outside the mask. Since the system of equations © is 
already over-determined, this additional information is read- 
ily extracted and reveals the temperature structure within 
the mask after the "reconstruction" . It should be noted that 
even a downgrade in the HEALPix resolution can lead to an 
information transfer by carrying out the downgrade before 
the mask is applied. 

The two masks used in this paper are the KQ85 (7yr), 
shown in figure QJ,, and the KQ75 (7yr) mask which are 
available a t the LAMBDA w ebsite. These are stored in the 
HEALPix JGorski et al.ll2005l ) format with a HEALPix reso- 
lution of iV s idc = 512. The reconstruction algorithm requires 
masks in lower resolutions of N s id« = 16 or 32, and a down- 
grade has to be carried out. A single downgraded pixel i 
contains iV tota i pixels of the higher resolution map and from 




(b) KQ85 (7yr) mask, iV aido = 16, x th = 0.5 
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Figure 1. The KQ85 mask of the WMAP 7 year data with a 
pixel resolution of N s ^ c = 512 is displayed in panel [l^,. The 
masked region is pictured in black. Figure [Tp shows the KQ85 
mask downgraded to a pixel resolution of -/V s ; dc = 16. All pixels 
with x(i) H 0.5, equation J5}, are considered as masked. 



these only iV nm pixel values are used in the averaging pro- 
cess, i. e. those that are not masked in the higher resolution 
map, in order to compute the value of the downgraded pixel. 
If the ratio 



x{i) = 



iV„ 



iV t , 



(5) 



is larger than a given mask threshold x t h £ (0.0,1.0), i.e. 
x(i) > Xth, the resulting pixel is not masked, and otherwise it 
is masked. The result of the KQ85 (7yr) nrask for AT s id e = 16 
and a; t h = 0.5 is shown in figureflp. It is obvious that the size 
of the masked domains in masks with N s id B < 512 depends 
on this mask threshold x t h- 

In the next sections we will compare results of CMB 
simulations of the ACDM concordance model with those of 
the ILC map. For this reason all maps of the ACDM con- 
cordance model are produced in a FWHM resolution of 1° 
and TVsidc = 512. But we also investigate these maps and the 
ILC map after an additional smoothing of e. g. 10° . In the 
following only the additional smoothing width is specified. 

Let us now return to the information transfer caused by 
the smoothing procedure and/or by carrying out the down- 
grade. The following simple numerical experiment presented 
in figures [2] and [3] reveals the information transfer. Figure[2^ 
displays the CMB simulation A of the ACDM concordance 
model after a smoothing of 10°. This map is downgraded 
from a pixel resolution of JV a id e = 512 to A/side = 16. Then 
the KQ85 mask in the pixel resolution N B id a = 16, shown 
in figure [TJa, is applied. For Z max = 10 the reconstruction al- 
gorithm ((3} is used to obtain the reconstructed map shown 
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(a) ACDM simulation A 



(a) ACDM simulation A 




-100 [JK 



(b) ACDM simulation B 



+ 100 [JK -10QmK. 



100 |JK 



(b) ACDM simulation B 




-100 [JK" 



+ 100 |jK -1 00 |JK- 



+100 mk 



(c) ACDM simulation A outside and B inside KQ85 mask 



(c) ACDM simulation A outside and B inside KQ85 mask 




-100|Jk' 



+ 100 [JK -100 mk 



+100 mk 



Figure 2. Figure|2]i shows the CMB simulation A of the ACDM 
concordance model with a pixel resolution of N s \^c = 512 and an 
additional smoothing of 10° . In figure[2p a second map (CMB sim- 
ulation B) using the same cosmological parameters is displayed. 
The map in figure[2j: results from the models used for figure[2k and 
[2p where the pixels outside the KQ85 mask are taken from simula- 
tion A and inside from simulation B at the resolution N B \^ e = 512. 
Thereafter a Gaussian smoothing of 10° is applied. 



Figure 3. The three reconstructed maps are displayed. The re- 
constructions are carried out for l max = 10. These maps result 
from the maps in figure [2^, b and c by downgrading these maps 
to a pixel resolution of A r s ; ( j c = 16 and thereafter applying the re- 
construction algorithm J3} to the data outside the mask shown in 
figure [TJ). Note that panel (a) and (c) should show the same sky 
map if the reconstruction algorithm would not use information 
from the masked region. However, there are significant difference 
especially near the Galactic plane. 



in figure [3K- Both the original and the reconstructed map 
agree within the mask reasonably well. The same procedure 
is repeated for a second simulation B where figure [2)3 shows 
the simulation and figure [3p the reconstruction. In the next 
step the pixels of simulation A are replaced at a pixel res- 
olution of A^idc = 512 within the KQ85 mask, shown in 
figure QJ,, by those of simulation B shown in figure [2]d. Af- 
ter this replacement the smoothing of 10° is applied. The 
last step transfers now the "wrong" information to the pix- 



els outside the mask. This smoothed map is shown in fig- 
ure [2J;. Downgrading this map to N a a e — 16 provides the 
data outside the mask which are used for the reconstruction. 
If the reconstruction would not use the information within 
the mask, the reconstructed map of figure [3^. should reap- 
pear. However, as revealed in figure [5J;, the reconstruction 
algorithm generates within the mask the main structures 
of simulation B, which is displayed in figure [2]d. This clearly 
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Figure 4. The leak of the ILC temperatures from regions within 
the mask to those outside due to the smoothing process. 



demonstrates the information transfer, so that one has to be 
careful in testing the reconstruction algorithm. This leads to 
the question whether the reconstruction can be carried out 
using only unsmoothed maps where no information about 
pixels within the mask is encoded outside. However, then 
stability difficulties arise as shown below. 

Another way to demonstrate the flow of information 
from regions within the mask to those outside the mask, is 
the following. At first the pixels of the ILC map outside the 
KQ85 mask are set to zero. Then a subsequent 10° smooth- 
ing shows how much information about the ILC pixel values 
inside the mask leaks to those regions outside. This leak of 
information is shown in figure [4] It is obvious that the main 
structures of the ILC map within the mask appear close to 
the boundary outside the mask. 

The paper is organised as follows. In section 2 the re- 
construction errors of the CMB temperatures are evaluated 
outside and inside the masks using different smoothings 
and different resolutions. The reconstruction method using 
the covariance matrix is compared with the direct inver- 
sion method. In section 3 the influence of the reconstruc- 
tion method onto the 2-point correlation function C(i?) of 
the CMB is investigated where the focus is on large scales. 
An integrated measure of C(#) serves in section 4 as a fur- 
ther tool to demonstrate the drawbacks of the reconstruction 
methods. Finally, in section 5 we summarise our results. 



2 STABILITY OF THE RECONSTRUCTION 

In the last section the quality of the reconstruction method 
is tested only visually. In order to quantify the accuracy 
of the method, we start with a temperature map <5Ttrue(i) 
which is considered as containing the true temperature in- 
formation also inside the mask. Here i denotes the index of a 
pixel within the map. The reconstruction algorithm ((3]) gets 
as input only the temperatures of those pixels lying outside 
the given mask and returns a map (5T rec (i,^max) for all pix- 
els. The resulting map depends on the multipole Z max up to 
which the reconstruction is carried out. The reconstruction 
accuracy can then be quantified by 



^Jrec {<> 



rcc lionaxj 
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Figure 5. The three upper curves show the normalised total 
error Q(/ m ax) of the reconstruction applied to the ILC (7yr) map 
using the KQ85 (7yr) mask in dependence on the multipole Z max . 
In addition to the original ILC map, the algorithm is applied to 
smoothed ILC maps with a smoothing of 300 and 600 arcmin at 
the resolution N s ^ a = 512. The reconstruction is carried out for 
the mask threshold £ t h = 0.5 using the resolution N a i^ c = 16. 
The three lower curves show the total error Q(£ max ) evaluated 
for the pixels outside the mask. 



where the prime means that the sum is restricted to the N m 
pixels within the mask. This gives the average error within 
the mask. The magnitude of S rec (I max ) depends on the nor- 
malisation of the map 5T tluc (i, l m&x ), which is obtained from 
the original map STtrue(i) by taking into account only the 
modes up to Z max . Note that the expansion ST tIU e{i, Li) is 
only possible for full sky maps since a spherical expansion is 
necessary. To normalise E rcc (/ max ) to a unit fluctuation, the 
reconstruction accuracy should be measured with respect to 
the mean temperature fluctuation 



^truc (,'m a x ) 



^-^V t 2 ruc (i,Z max ) , (7) 



where the sum with the two primes is restricted to the N 
pixels outside the mask. The normalised total reconstruction 
error is then given by the quotient 

-2-Tec ^'m a x ) 



Q(ln 



CT tr 



s 



(8) 



rue \ t-niax^ 



(6) 



A value of Q = 1 means that the error is as large as a typical 
temperature value, i. e. the reconstruction is useless. 

We apply the normalised accuracy measure (JS]) to the 
ILC (7yr) map and to simulations of the ACDM concordance 
model where the cosmological parameters of the WMAP 
team are use d (Tab le 8, column "WMAP+BAO+J7 " in 
Ijarosik et al.l (|2010j)). These maps with a resolution of 
-Wsidc = 512 are smoothed to a Gaussian width 6 = 1°, 
and in the following only a smoothing in addition to this 
8 = 1° smoothing is given. After the smoothing the maps 
are downgraded to lower resolutions of A^idc = 16 or 32. 
Then the downgraded masks are applied in this lower res- 
olution, which gives the maps which serve as the input for 
the reconstruction algorithm. 
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(a) inside KQ85 mask 
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(b) outside KQ85 mask 



+ 1.00 
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Figure 6. The upper panel shows the local reconstruction error 
J9jl of the ILC (7yr) map within the KQ85 mask corresponding 
to the total reconstruction error Q(10) = 0.62 of the ILC map 
as can be read off from figure [5] at Z max = 10. The colour bar is 
truncated such that q-values less than -1.0 are pictured as dark 
blue and larger than 1.0 as dark red. In the lower panel the same 
quantity is pictured outside the mask. Here typical g-values lie 
between -0.2 and +0.2 as can be read off in figure \5\ For this 
reason the maximum and minimum values are truncated at 0.2 
and -0.2, respectively. 



In figure [5] the three upper curves reveal the total re- 
construction error (|8]) for ILC (7yr) map using the KQ85 
(7yr) mask for different smoothings. As outlined above the 
algorithm depends on the multipole Z max up to which the 
multipoles 1 = 0,..., / max are to be reconstructed. The ac- 
curacy measure (JSJ) refers to the maps expanded only up to 
I = Imai. Even for very low values of Z max < 10 the total er- 
ror Q(£ max ) is above 0.4, i. e. the temperatures reconstructed 
within the mask differ from the true ones by one half of the 
mean temperature fluctuations. This demonstrates the ac- 
curacy problems of the algorithm. It is seen that the error 
decreases with increasing Gaussian smoothing because the 
information leak is getting stronger with increasing smooth- 
ing. The three lower curves shown in figure [5] are the total 
errors outside the mask. These are computed using equation 
([6]) but with the modification that the primed sum runs only 
over those pixels i outside the mask. It is seen that with in- 
creasing value of £ max the accuracy gets better although the 
error is surprisingly large for small reconstruction multipoles 
' max - With additional smoothing the total reconstruction er- 
ror of the ILC map has the tendency to fall, but without ad- 
ditional smoothing the total error is nearly constant about 
0.2. 



ILC. Ojirernm 

ILC, 300 arcmiti 

ILC, 600 STcmin 




Figure 7. The same analysis for the normalised total error 
Q(Z max ) as in figure \5\ is presented, but now the larger KQ75 
(7yr) mask is used instead of the KQ85 (7yr) mask. 



Up to now only the global accuracy measure Q is dis- 
cussed. But in order to find the regions which are most prone 
to reconstruction errors, we now introduce the reconstruc- 
tion accuracy of single pixels which is obtained by taking 
the difference between ST TCC (i, Z max ) and <5T tr uo(i, Z max ). Nor- 
malising this difference with o"true(^m ax ) results in the local 
reconstruction error 



q(i,l n 



01 rcc \1, 'n 



5T trU e{i,L 



Ctr 



.(J 



(9) 



rue I "ma.x.1 



To expose the regions with the largest errors, the local 
reconstruction error ([9]) of the ILC map (without additional 
smoothing) is pictured in figure [6] for Z max = 10. The errors 
within the mask are shown in figure [6^. In this figure dark 
areas display errors with |g| ^ 1. Such values of q occur 
near the galactic centre where the mask is very extended. 
But errors larger than 1 also occur antipodal to the galactic 
centre. The errors outside the mask are shown in figure [6p. 
Based on figure [5] one would expect that a typical value of 
the local reconstruction error @ is about ±0.2. The colour 
band is truncated at those values of |g| such that values of 
\q\ ^ 0.2 are pictured as dark blue or dark red. It is seen 
that large errors outside the mask appear near the masked 
region. A further observation is that with increasing Z max the 
most significant errors approach the masked region. Then 
the values of the maximum errors increase, although the 
areas with \q\ ^ 0.2 are getting smaller. 

Since the reconstruction error depends on the size of 
the mask, it is interesting to consider the larger KQ75 (7yr) 
mask which again is applied to the ILC (7yr) map. In figure 
[7] the three upper curves display the corresponding total re- 
construction errors ([8)l within the mask. Even for very low 
values of Z max < 5 the total error Q(Z max ) is above 0.5, i. e. 
the temperatures reconstructed within the mask differ from 
the original ones by more than one half of the mean tem- 
perature fluctuation. For Z max ^ 7 the total error Q(7 max ) 
is even above 1.0. It is seen that the error decreases with 
increasing Gaussian smoothing because of the information 
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(a) inside KQ75 mask 
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-1.00 



(b) outside KQ75 mask 



+ 1.00 
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Figure 8. The upper panel depicts the local reconstruction error 
l|9)l of the ILC (7yr) map within the KQ75 mask corresponding 
to the total reconstruction error Q(10) = 1.86 as figure [7] reveals 
at /max = 10. In the lower panel the same quantity is pictured 
outside the mask. The colour bands are truncated in the same 
way as in figure [6] i.e. at ±1.0 and ±0.2 in the upper and lower 
panel, respectively. 



leak, but nevertheless the total error Q(/ ma x) is large for a 
smoothing with 600 arcmin. The three lower curves shown 
in figure [7] are the total errors outside the mask. It is seen 
that with increasing value of £ max the accuracy improves 
although the error is surprisingly large for small reconstruc- 
tion multipoles Z max . With additional smoothing the total 
reconstruction error of the ILC map has the tendency to 
fall, but without additional smoothing the error is nearly 
constant at about 0.25. 

To reveal the locations of the most significant errors in 
the case of the KQ75 (7yr) mask, the local reconstruction 
error of the ILC map (without additional smoothing) is pic- 
tured in figure [H] for Z max = 10. The errors within the mask 
are shown in figure [8^. The dark areas display errors with 
|g| JS 1 which are much more widespread than for the KQ85 
(7yr) mask. The errors outside the mask are shown in figure 
[8p. To locate large errors we have pictured values of |q| > 0.2 
as dark blue or dark red. The errors outside the mask are 
only slightly larger than in the case of the KQ85 (7yr) mask. 

Compared to the KQ85 (7yr) mask, the accuracy prob- 
lems of the algorithm applied to the larger KQ75 (7yr) mask 
are much stronger, and we conclude that this mask is too 
large in order to allow a reconstruction for Z max ^ 7, and 
even the results for Z max < 7 should be considered critical. 
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Figure 9. In the upper panel the solid curve shows the normalised 
total error Q(£ ma x) of the reconstruction inside the KQ85 (7yr) 
mask applied to the ILC (7yr) map smoothed with 600 arcmin 
in dependence on the multipole Z max . The reconstruction is car- 
ried out for the mask threshold x t \^ = 0.5 using the resolution 
^sidc = 16. The mean value of the total reconstruction error 
calculated from 1000 CMB simulations is displayed as a dashed 
curve. The range containing 66.6 % of these models is pictured as 
dark grey band. The distribution of all errors of these 1000 CMB 
simulations is given as a light grey band. The lower panel shows 
the corresponding total reconstruction errors outside the mask. 



For this reason we concentrate us in the following on recon- 
structions using the smaller KQ85 (7yr) mask. 

Now we compare the total reconstruction error of the 
ILC map with those of 1000 CMB realizations of the ACDM 
concordance model. In figure [5] this comparison is shown 
with the additional smoothing of 600 arcmin and in figure 
1101 without smoothing. A resolution of iV a idc = 16 and a 
mask threshold x t h = 0.5 is used. As can be seen in the 
upper panel of figure [9] the total reconstruction error of the 
ILC map within the mask is located for Z max ^ 8 outside 
the dark grey band containing 66.6% of 1000 models, but 
between the maximum and minimum errors of these mod- 
els which domain is depicted as the light grey band. For 
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I mal > 8 the total reconstruction error of the ILC map is 
rather typical. It should be noted that the computation of 
the ILC reconstruction error assumes that the ILC map con- 
tains the true values inside the mask which, of course, needs 
not to be the case. But as the comparison with the light grey 
band shows, the total reconstruction error is always smaller 
than the extreme errors that occur among the 1000 simula- 
tions. The total reconstruction error from these 1000 models 
with an additional smoothing of 600 arcmin reveals that the 
reconstruction of the CMB within the mask is impossible for 
I max > 16. In figure [9p the total reconstruction error of the 
ILC map outside the mask is compared with those of the 
1000 ACDM models. For 2 max ^ 8 the total reconstruction 
error of the ILC map is again outside the dark grey band 
containing 66.6% of the simulations. Furthermore, with in- 
creasing £ max the total reconstruction error outside the mask 
vanishes if the additional smoothing of 600 arcmin is applied. 
Figure [10] displays the total reconstruction error with- 
out the additional smoothing of figure [9] The total recon- 
struction error within the mask of the ILC map and of 
the 1000 simulations is larger without additional smoothing. 
With increasing value of Z max the reconstruction of the tem- 
perature fluctuations within the mask deteriorates increas- 
ingly, such that the reconstruction is doubtful for / max > 13. 
The mean values of the total reconstruction error outside 
the mask of the 1000 ACDM simulations are now nearly in- 
dependent from / max at a value about 0.15. The width of the 
distribution of the total reconstruction error decreases with 

'm aX - 

Figures [6] and [8] depict the local reconstruction error 
q{i, £ max ) for the ILC map, where one has to assume that 
the ILC map contains the true pixel values within the mask. 
The locations of extreme values of q(i, i max ) without this 
assumption are obtained from the 1000 ACDM simulations. 
The standard deviation of q(i, £ max ) calculated from these 
1000 models (without additional smoothing) is pictured in 
figure im for £ max = 10. The upper panel shows the standard 
deviation within the mask and the lower one outside the 
mask. Within the mask the standard deviation is extremely 
large near the centre of the galaxy where a huge area is 
masked. In general, one finds that the larger the masked 
domain is, the larger is the standard deviation of q(i, Z max ) in 
the centre of such a domain, and accordingly, the larger is the 
expected local reconstruction error for an individual map. 
Outside the mask, extreme standard deviations of q(i, Z max ) 
occur near the boundary of the mask. The comparison of 
the local reconstruction error in figure [6] with this standard 
deviation shows that antipodal to the centre of the galaxy 
the local reconstruction error is with about 4rr untypically 
large. This could be a hint for residual foreground in the 
ILC map. 

The reconstruction error depends also on the resolution 
of the maps. In order to address this point, figure [12] com- 
pares the total reconstruction errors Q(7 max ) obtained from 
maps with HEALPix resolutions Aside = 16 and 32. The 
distribution of Q(/ max ) of 1000 CMB simulations as well as 
the curves belonging to the ILC maps are shown. The maps 
are smoothed to 600 arcmin and the KQ85 (7yr) mask with 
a mask threshold x t h = 0.5 is used. One observes that the 
reconstruction performs at the higher resolution A s id c = 32 
a bit better than at the lower resolution Aside = 16. This be- 
haviour takes place within and outside the mask in a similar 
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Figure 10. The same analysis of the reconstruction error Q(Z max ) 
as in figure [9] is shown, but now without smoothing. 



way. As previously discussed, it is important to carry out the 
reconstruction without a smoothing at intermediate steps. 
For that reason, figure [TH1 shows the total reconstruction er- 
rors Q(Z max ) without additional smoothing. Aside from this 
difference the analysis is the same as in figure [j"2l Omitting 
a smoothing procedure leads to an overall increase in the to- 
tal error Q(7 max ). In addition, without smoothing the recon- 
struction at the higher resolution N a ide = 32 performs only 
marginally better than at the lower resolution A^ide = 16. 
Thus, reconstructions of CMB maps at higher resolutions 
A s idc do not yield a crucial inrprovement if a smoothing has 
to be avoided. 

Now we will discuss the influence of the mask thresh- 
old Zth defined in equation ([5}. The lower the value of Xth, 
the more pixels are excluded in the neighbourhood of the 
boundary of a masked domain. The total reconstruction er- 
rors Q(7 max ) of 1000 CMB simulations of the ACDM con- 
cordance model for x t h = 0.1, 0.5 and 0.9 at a resolution 
A s idc = 16 (without additional smoothing) are compared 
in figure [T4l For xth = 0.1, 0.5 and 0.9 the reconstruction 
method uses 84.9%, 80.8% and 70.4% of the pixels of a map, 
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Figure 11. The standard deviation of the local reconstruc- 
tion error J9} is displayed for /max = 10 calculated from 1000 
ACDM maps. The upper panel depicts the standard deviation of 
q(i, /max) within the mask, where values 9^1 are pictured white. 
In the centre of the galaxy the standard deviation has values as 
large as 2.12. At the lower panel the same quantity is pictured 
outside the mask, but now the maximum value is truncated at 
0.2. Pixel values larger than 0.4 occur. 



respectively. Thus the KQ85 (7yr) mask for x t ^ — 0.9 at a 
resolution -Nside = 16 excludes a region almost as large as 
the initial KQ75 (7yr) mask, and for x t h = 0.5 the excluded 
region is even smaller than the initial KQ85 (7yr) mask. 
Since the reconstruction accuracy decreases by increasing 
the mask and because a larger value of a?th increases the 
mask, the reconstruction has more difficulties for larger Xth 
as can be seen in figure [T4l One observes in the upper panel 
that the differences in the total reconstruction errors within 
the mask for the three thresholds x%h are relatively small as 
long as Z max < 5. But with increasing multipoles Z max the 
total reconstruction error increases strongly with i t h- The 
total reconstruction error within the mask for x t ^ = 0.1, 
0.5 and 0.9 is typically larger than 1 for Z ma x > 15, 13 and 
11, respectively. These represents the largest multipoles Z ma x 
for which a reconstruction can be carried out in order to 
avoid constructing random patterns. But for cosmological 
parameter extractions, even values of Q(7 m ax) near one can 
be insufficient, since this implies an total error of a typical 
temperature fluctuation. 

Up to now, all reconstructions have used the covariance 
matrix (J4} in equation Q. Thus, we finally compare the 
reconstruction results by using the covariance matrix with 
those of the method of direct inversion, where A is the unit 
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Figure 12. The normalised total error Q(Z m ax) is computed for 
the two different resolutions of A s id e = 16 and A s id e = 32. The 
total error Q(Z max ) is obtained from the ILC (7yr) map together 
with KQ85 (7yr) mask using a threshold x t ^ = 0.5. A smoothing 
of 600 arcmin is applied. The total error for the ILC (7yr) map 
for Agije = 16 is plotted as a solid curve and for A s id e = 32 as a 
dash-dotted curve. This error is compared with the corresponding 
mean error obtained from 1000 ACDM simulations which is dis- 
played as a dashed curve and as a dotted line, respectively. The 
distribution of the errors of the 1000 CMB simulations is given 
as a light grey band and as a shaded band for A s id e = 16 and 
A s ido = 32, respectively. The upper and lower panels show the 
total reconstruction errors inside and outside the mask. 



matrix. Choosing the unit matrix renders the method inde- 
pendent of an underlying cosmological model, see equation 
(J4|. In figure [151 this comparison is carried out for the 1000 
ACDM simulations using the KQ85 (7yr) mask for a mask 
threshold x t h = 0.5 at a resolution of iVside = 16. The re- 
construction using the covariance matrix works better than 
the method of direct inversion for ( m « < 8. For Z max > 9 the 
difference between the two methods is marginal. Thus, for 
the reconstruction of the lowest multipoles the algorithm 
using the covariance matrix is preferable, whereas for the 
larger multipoles the direct inversion can be chosen due to 
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Figure 13. The normalised total error Q(( m ax) is computed from 
the same data as in figure [12] The sole difference is that no addi- 
tional smoothing is applied here, which leads to larger values of 
Q('max) compared to figure fl2l Furthermore, a higher resolution 
leads only to a marginal improvement. 



its smaller computational effort. In the case of the method of 
direct inversion, the mean value of the total reconstruction 
errors of 1000 simulations at a resolution of iV S ide = 256 
are also calculated, which does not result in a crucial im- 
provement compared to the reconstruction at a resolution 
of JVsidc = 16. The similarity of both variants of the recon- 
struction algorithm is due to the fact that the covariance 
matrix is already nearly diagonal. 

The reconstructions of the 1000 ACDM simulations and 
of the ILC map are compared with the multipole expansions 
of the full initial maps up to the reconstruction value of Z ma x- 
Thus, this comparison stringently needs the full sky map in 
order to carry out the multipole expansion. In the case of 
the simulations we know that the full initial maps contain 
the so-called true CMB temperature fluctuations, and the 
expansion can be done reliably. The case of the ILC map 
is more involved since the ILC map contains residual noise 
and foreground. In particular, within and near the mask one 
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Figure 14. The mean value of the total reconstruction errors 
Q(^max) calculated from 1000 ACDM simulations using the KQ85 
(7yr) mask for the mask thresholds x^ n = 0.1, 0.5 and 0.9 at a 
resolution N s ^ e = 16 are compared. In the case of x t ^ = 0.9 the 
distribution of all reconstruction errors is given as a grey band. 
The upper and the lower panel show the total reconstruction er- 
rors computed within and outside the masks, respectively. 



would expect remaining foreground to some extend. The so- 
called true temperature fluctuations of the ILC map are not 
necessarily the genuine CMB temperature fluctuations. One 
has to keep this fact in mind, when one tries to weight the 
quality of the ILC reconstruction on the basis of the total 
reconstruction errors Q(Z ma x) of the simulations. Thus it is 
possible that the reconstruction of the CMB of our Universe 
from the ILC map is even worse than the result in figure [TOl 
suggests. 



3 CORRELATION FUNCTION C(0) 

In this section we address the crucial question whether the 
sky reconstructions can reliably be used for the estimation of 
cosmological statistics. The correlation function C(i?) is the 
suitable statistic with respect to large angular scales # which 
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Figure 15. The total errors Q(£ max ) of the reconstruction using 
the covariance matrix Q (cov) are compared with the total errors 
of the direct inversion (di). Both methods use the KQ85 (7yr) 
mask with a mask threshold a; t h = 0.5 at a resolution N s ^ c = 16. 
The direct inversion is also computed for -/V^e = 256. The total 
error distribution of the 1000 CMB simulations is given as a dark 
grey band for the algorithm using the covariance matrix and as 
a light grey band for the direct inversion (N a ^ e = 16). In the 
upper panel the total reconstruction errors Q(£ max ) within the 
mask are given, and it is seen that for £ max ^ 8 the covariance 
method possesses a smaller error distribution. The lower panel 
shows the corresponding total reconstruction errors outside the 
mask. 



are relevant for the reconstruction algorithm. This leads to 
the question of how sensitive the correlation function C(#), 
equation {TJ, varies with respect to the different sky recon- 
structions. At first one has to determine the minimal value of 
the reconstruction parameter / max which is needed for the 
reconstructed map in order to give a correlation function 
C($) in agreement with the one obtained from the initial 
map with full resolution. Therefore, the 2-point correlation 
functions computed from the multipoles up to 2 max = 5, 10 
and 15 are compared with the correlation function obtained 
from the full ILC map with maximal resolution in figure [16] 



Figure 16. The correlation functions of the full ILC map calcu- 
lated from the multipoles up to Z max = 5, 10 and 15 are compared 
with the correlation functions of the full ILC map with maximal 
resolution (max). 



From this figure one can read off that at least a reconstruc- 
tion parameter / max > 10 is needed in order to have a good 
estimation of the correlation function at large scales. The 
investigations of the previous section lead to the conclusion 
that the KQ75 mask is already too large in order to allow 
a reconstruction for such large values of Z max . For this rea- 
son the following computations are based all on the KQ85 
mask with a mask threshold a;th = 0.5 at a resolution of 
-Wsidc = 16. Furthermore, the reconstruction algorithm us- 
ing the covariance matrix is applied in the following, since 
it gives for Z max > 10 almost the same results as the direct 
inversion method. 

In figure [17] three correlation functions C(l?) are shown 
which are calculated from the reconstructed ILC (7yr) map, 
from the full ILC map and from the data outside the KQ85 
mask. The three panels show the results for the reconstruc- 
tion parameters Z max — 10, 12 and 15. Here the ILC (7yr) 
map with an additional smoothing of 600 arcmin is used. 
The reconstructed correlation function is given with ler and 
2cr deviations of the reconstruction errors which are esti- 
mated by 



(t(0,L 






C tluc ($,lmax)) 



N 



(10) 
where the index k counts the N = 1000 CMB simulations 
of the ACDM model. This quantity takes into account only 
the error of the reconstruction of pure CMB maps. Addi- 
tional contributions due to residual foreground and detector 
noise in the ILC map propagate through the reconstruction 
algorithm and are not considered here. In all three cases 
(^ max = 10, 12 and 15) the correlation function of the recon- 
structed ILC map agrees well within the 2a errors with the 
correlation function of the full ILC map and less with the 
correlation functio n computed outside the mask. This ob- 
servation is used bv lEfstathiou et al.l (120091 ) as an argument 
in favour of the statement that the "true" correlation func- 
tion is the one obtained from the full ILC map. However, as 
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already emphasised, the large additional smoothing of 600 
arcmin transfers information from the masked region into 
the data used by the algorithm. For this reason this result 
is not surprising, and the reconstructed correlation function 
is very questionable since it is based on the pixels which 
should have been omitted. 



To avoid the problem of information transfer by 
smoothing the correlation function of the ILC (7yr) map 
is also investigated without additional smoothing. In figure 
[T5]the corresponding correlation functions are plotted. The 
error (|10|) is significantly larger now. Without additional 
smoothing the differences between the correlation function 
of the reconstructed and of the full ILC (7yr) map are larger, 
but both agree within the 2<r errors. This could lead to the 
conclusion that the reconstruction works without additional 
smoothing and without a corresponding information trans- 
fer, but for / m ax = 10 — 12 the errors are large compared to 
the case with additional smoothing as it is shown in figure 
1171 For / max > 15 the errors of the reconstructed correlation 
function are too large to allow any conclusion. Due to these 
large errors the correlation function of the reconstructed ILC 
map also agrees on the same level with the correlation func- 
tion obtained solely from the data outside the mask. Thus, 
one cannot differentiate between these two cases. The re- 
constructed correlation function shows sometimes a better 
agreement with the correlation function resulting from the 
masked ILC map and sometimes with the correlation func- 
tion of the full ILC map. The quality of this agreement de- 
pends on $ and Z max . Thus, the reconstructed correlation 
function favours neither the one of the masked nor the one 
of the full ILC map. Due to the large errors the reconstructed 
correlation function is uncertain by at least 100/iK 2 and thus 
unsuited for the comparison with cosmological models. 
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The above discussion puts forward arguments against 
the reconstruction method and favours methods which use 
only the data outside a given mask. To provide a firm footing 
for the latter, the influence of the KQ85 and KQ75 masks 
onto the ensemble average and the cosmic variance of the 
correlation function C($) is now investigated with respect 
to the ACDM model. Both quantities are shown in figure fUJl 
for an ensemble of 100 000 CMB simulations of the ACDM 
model. The ensemble average and the standard deviation 
is computed from the correlation functions C(i?) obtained 
from the data of full maps, outside the KQ85 mask and 
outside the KQ75 mask. A resolution of A^ide = 128 and a 
mask threshold Kth = 0.5 is used. The ensemble averages are 
identical in all three cases. The standard deviation slightly 
increases with the size of the mask. The smallest 1<t standard 
deviation is obtained by using no mask at all, the next larger 
deviation belongs to the KQ85 mask whereas the largest de- 
viation is due to the larger KQ75 mask. The increase of the 
standard deviation for the masked data is, however, small 
compared to the uncertainty resulting from the reconstruc- 
tion method for the same mask. This is a further argument 
for having more confidence in the correlation function C($) 
computed from the data outside the mask than in the one 
obtained from a reconstructed full map. 
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Figure 17. In panel (a), (b) and (c) the correlation function 
C(i9) of the reconstructed ILC map (cov) is compared with C(i9) 
of the full ILC map (full) and C(i9) computed outside the mask 
(mask) for ( max = 10, 12 and 15, respectively. Here the ILC map 
is additionally smoothed to 600 arcmin. The reconstruction is 
based on the KQ85 mask (x t ^ = 0.5 and N s ^ e = 16). The dark 
grey and light grey band show the la and 2a band of a(i), /max), 
eq. llOt . calculated from 1000 ACDM simulations, respectively. 



© 2010 RAS, MNRAS OOO.rTHTsl 



12 R. Aurich and S. Lustig 



C(«) 
ax 

200 



-400 



ILC, arcmiu, fuJ] 
TLC, Oarcman, cov 
[LC, fljcmiii, mflsfc 




100 000 ACDM. full 

I DO COO ACDM, KQ85 mask 

100 000 ACDM, KQ75 monk 




180 ■& 



400 

[ M K 2 ] 



V.I i 



-400 



- LLC, aiuuiiij , full 
• ILC, arcmin | «jv 
ILC I ' urciuin, musk 




300 

'I 

•sec 

-400 


111 

I'll 1 

llll \ 

1 : , 11 '. 

111 \ \ 
\i\ \ \ 

U\\\ /\ 


. | i i | i . 




ILC, aremja, full 

ILC, arcmin,. cov 


/" 


/ 

/ 




\X/ ^ 

/ 

_ i'A^J 'max == J-0 

/ 

. /. 1 . . 1 . . 1 , 


//~X \ \ > 

•> / \ \ ■■ 

/ \ \ \ '■-- 

\ \ 



U.I 



•III 



:;:i 



\>ii 



I.",. I 



.SI I 



1? 



Figure 18. The same computations are displayed as in figure [TTl 
but here the ILC (7yr) map is used without additional smoothing 
which leads to significantly larger errors. 



Figure 19. The mean value of the correlation function C(#) of 
full maps computed from 100 000 CMB simulations of the ACDM 
model is compared with the mean values of C(i?) using the KQ85 
(7yr) mask and the KQ75 (7yr) mask. The mean values cannot 
be distinguished in this plot. All calculations are carried out with 
a resolution of N s -^ c = 128 and a mask threshold of x t ^ = 0.5. 
The light grey band gives the la standard deviation calculated 
from full maps of the 100 000 CMB simulations. The medium grey 
and the dark grey band correspond to the la standard deviations 
calculated from the same 100 000 maps, but now outside the KQ85 
and outside the KQ75 mask, respectively. 



4 THE LARGE ANGULAR SCALE POWER 

The main difference between the ILC correlation functions 
C(#) computed from the data outside the KQ85 mask and 
from all pixels is that the former displays much less power 
at angular scales above 60° than the latter. To describe this 
observation the S^O ) statistic 



S^O ) := 



cos(60°) 



dcostf |C(tf)| 



(11) 



has been introduced (jSpergel et al.l 120051 ) . which measures 
the power on scales larger than 60°. In the following the 
5(60° ) statistic is discussed as a function of the reconstruc- 
tion parameter £ max - At first the results obtained from maps 
with an additional smoothing of 600 arcmin are presented. 
In figure [20k the 5(60°) statistic is computed from the full 
ILC map, from the data outside the mask and from the re- 
constructed ILC map. It is striking to see that in the case 
of the reconstructed map the S(60°) statistic does not con- 
verge to a stable value as it is the case for the full ILC map 
as well as for the masked ILC map for <! max > 6. This em- 
phasises the uncertainty of reconstructed maps despite the 
information transfer by smoothing the map. In figure 120b 
the mean value of the S(60°) statistic, calculated from 1000 
CMB simulations of the ACDM model, is displayed again 
with an additional smoothing of 600 arcmin. As in figure 
120b . the statistic is plotted for full maps, for reconstructed 
maps and for masked maps. Since the ACDM model has 
significantly more power at scales above ■& — 60°, the scale 
of the ordinate now differs from that of figure 120b . This 
should be noted by comparing both panels. For Z max < 10 
the mean value of the S^O ) statistic of the reconstructed 
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Figure 20. In panel (a) the S(60°) statistic of the reconstructed 
ILC (7yr) map (cov), the full ILC (7yr) map (full) and the ILC 
(7yr) map outside the applied mask (mask) is plotted as a func- 
tion of /max- Here the ILC map is additionally smoothed with 
600 arcmin. The reconstruction is carried out applying the KQ85 
(7yr) mask for the mask threshold x tn = 0.5 using the resolution 
^side = 16. In panel (b) the mean values of the 5(60°) statistic 
calculated from 1000 CMB simulations of the ACDM model are 
plotted. 



maps differs only marginally from that of the full maps, but 
for Z max > 10 the differences are pronounced. If the recon- 
struction method would work perfectly, the mean values of 
the 5(60°) statistic of reconstructed maps and of full maps 
should be identical, but there is a systematic overestimation 
of power in the reconstructed maps. This clearly hints to a 
systematic error by using the reconstruction method. The 
overestimation of the mean value of the ,5(60° ) statistic ob- 
tained from 1000 simulations again stresses the problems of 
the reconstruction method for Z max > 10. 

Furthermore, figure [20b reveals that the mean value of 
the 5(60° ) statistic increases by 10% by using only the data 
outside the KQ85 mask compared to the full pixel input. It 
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Figure 21. The same quantities as in figure l20l are displayed for 
the ILC (7yr) map in panel (a) and for 1000 ACDM simulations 
in panel (b), but now they are calculated from the maps without 
additional smoothing. 



is interesting to note that the reverse behaviour occurs in the 
case of the ILC map. There, the 5(60°) statistic obtained 
from the masked map is below that of the full map. Using 
the above observation based on 1000 simulations, one would 
rather expect that the 5(60°) statistic of the full ILC map 
should lie 10% below the one of the masked map. The ILC 
map displays in this respect an untypical behaviour. 

In figure [21] the same quantities are displayed as in fig- 
ure 1201 but now calculated without additional smoothing 
of the maps. Compared to the previous case with smooth- 
ing, the 5(60°) statistic of the reconstructed ILC map fluc- 
tuates even stronger for J max > 6 which emphasises the 
strong uncertainty of the reconstruction without the infor- 
mation transfer caused by smoothing. For Z max > 15 the 
5(60°) statistic of the reconstructed ILC (7yr) map yields 
extremely large values which corresponds to the overestima- 
tion of the power of the correlation function at large scales. 
The mean value of the 5(60° ) statistic calculated from 1000 
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CMB simulations of the ACDM model is plotted in figure 
121b . Again the overestimation of power in the case of re- 
constructed maps is much more pronounced as in the case 
without an additional smoothing. Now this overestimation 
increases very fast already for Z max > 12 which shows that 
for these values of l max the reconstructed results definitely 
contain significant systematic errors. For individual maps 
the errors can be significant for even lower values of Z max as 
shown in figure |2"TR for the ILC map. 

The figures [20l and [2"T1 also reveal the fact that the low 
value of S(60°) of the ILC map is rather unusual. As can be 
read off from the figures this value is much smaller than the 
mean value due to the 1000 ACDM simulations. The 5(60° ) 
statistic of the ILC map is now compared with those of 
100 000 ACDM simulations in order to find the frequency of 
simulations with even lower power than contained in the ILC 
map, i. e. which satisfy the condition (jSpergel et al.l 120031 ; 
ICopi et al.ll2010r) 

Sacdm(60°) < Silo (60°) . 

This condition is satisfied by only 5244, 54 and 16 models 
out of the 100 000 CMB simulations where S^O ) is com- 
puted from data of the full map, outside the KQ85 mask, and 
outside the KQ75 mask, respectively. As discussed above the 
values of 5(60° ) of the ILC map is lower for the masked map 
than for the full map which is unusual. How unusual this 
behaviour is, becomes evident by the fact that none of the 
100 000 ACDM simulations simultaneously satisfy the condi- 
tions Sacdm(60°) > Silo (60°) for data of the full map and 
5acdm(60°) < Silc(60°) for data outside the KQ85 mask. 
Using instead of the KQ85 mask the KQ75 mask leads to 
only 2 models out of 100 000 CMB simulations which sat- 
isfy both conditions simultaneously. This is a clear hint that 
the ILC map shows an untypical behaviour compared to the 
ACDM concordance model. 



5 SUMMARY 

In this paper we have investigated the reconstruction of pixel 
values within masks using inversion methods. We have com- 
pared the algorithm using the covariance matrix with the 
direct inversion method. The comparison shows that for the 
reconstruction of multipoles £ max < 8 the algorithm using 
the covariance matrix gives slightly smaller errors than the 
direct inversion method. For this reason we focus on the 
reconstruction method which uses the covariance matrix. 

An important point of our analysis is tha t the addi- 
tional smoothing of 10° of the ILC map in Efstathiou et al. 
(2009) transfers information from pixels inside the mask to 
pixels outside. This should not happen since the pixel val- 
ues inside the mask are considered as insecure and should be 
ignored. The smoothing, however, leads to a quantitatively 
better reconstruction of the ILC map, if the masked domain 
is considered as containing the true pixel values which have 
to be recovered by the reconstruction. But the contribution 
of information from masked pixels makes a reconstructed 
map obtained from data with additional smoothing unac- 
ceptable in a realistic application. 

The errors of the pixel values of the reconstructed 
map are compared with the mean temperature fluctuation 
0"truc(7max) defined in equation ©. In the case of the KQ75 



(7yr) mask the reconstructed pixels inside the mask have 
errors larger than fftme(imax) already for i max > 6 such that 
the reconstruction is unusable as revealed by figure [7] Even 
for smaller values of Z max the errors are at least of the order 
of 5<7truc(lmax). The culprit for this negative result is the size 
of the KQ75 (7yr) mask being too large for a reconstruction. 

In the case of the smaller KQ85 (7yr) mask the recon- 
struction without additional smoothing leads for Z max > 13 
to errors larger than o" trU e(^ma X ) such that a rough recon- 
struction can be done for Z max around 10. However, even 
with this restriction, errors larger than ^cVueCmax) occur, 
see figure [5] This might be too large for cosmological appli- 
cations. It is also checked that refining the resolution from 
Aside = 16, which is used in the above computations, to 
A a ido = 32 does not improve the reconstruction. 

Furthermore, the behaviour of the temperature 2-point 
correlation function C($) is analysed with respect to the var- 
ious reconstructed maps, i. e. with and without additional 
smoothing, for the ILC (7yr) map and for CMB simulations 
of the ACDM concordance model. The correlation function 
C(i9) is computed from the full ILC map, from the data out- 
side the KQ85 mask, and from the reconstructed ILC maps 
which depend on Z max . In order to get a good approximation 
of the correlation function C(i9) on large angular scales, the 
multipoles up to Z max ~ 10 are needed. The application of 
the large KQ75 mask is therefore excluded. The errors of 
the reconstructed 2-point correlation function are estimated 
by using 1000 CMB simulations of the ACDM model. Since 
this error estimation is based on simulations which contain 
only the pure CMB signal, no further errors are taken into 
account e. g. resulting from residual foreground and detec- 
tor noise. Therefore, the error estimation of the correlation 
function of the ILC map gives only a lower bound of the 
errors expected in a genuine application. 

As discussed above the reconstructed ILC map with ad- 
ditional smoothing of 600 arcmin uses information from in- 
side the mask. This information transfer leads in turn to 
underestimated errors in the 2-point correlation function. 
Furthermore, the similarity of C(i?) computed from the full 
ILC map and from the reconstructed map is due to this 
information transfer and not the achievement of the recon- 
struction algorithm. 

Without additional smoothing the ILC correlation func- 
tion C(t9), computed from the reconstructed map, does not 
match the one of the full map better than the one of the 
masked map within the estimated errors. This is caused by 
the large errors which increase with Z max . For 10 < Z max < 12 
the reconstructed 2-point correlation function is uncertain 
by 100/iK 2 and for even larger values of Z max useless as re- 
vealed by figure [TBI 

This uncertainty is also reflected in the 5'(60°) statistic 
which integrates the power of the 2-point correlation func- 
tion on angular scales larger than 60°. The S'(60°) statistic 
of the reconstructed ILC map does not converge to a sta- 
ble value for Z max > 6 as it is the case for the full and for 
the masked ILC map as displayed in figure l2"Tk . In addition, 
the mean value of the 5(60° ) statistic, calculated from 1000 
reconstructed CMB sky maps, reveals the instability of the 
algorithm for Z max > 12, where the power strongly increases, 
see the dashed curve in figure 121b . This demonstrates the 
inability to reconstruct the temperature correlations. 

The conclusion of this paper is that by using a realis- 
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tic mask for the WMAP data the reconstruction algorithm 
([3]) does not work well enough to obtain a prediction of the 
2-point correlation function. For this reason a cosmological 
analysis of the WMAP data should use only data outside a 
mask, eg. the KQ85 (7yr) or KQ75 (7yr) mask. The Planck 
satellite measures the sky at more different frequencies than 
the WMAP satellite. This should allow a more secure re- 
duction of the foreground in the CMB map. It is expected 
that the masked region in the Planck data is significantly 
smaller. It might be small enough in order to allow an ac- 
ceptable reconstruction of the CMB within the mask which 
could provide a 2-point correlation function useable at large 
scales. 
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